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Modulation of the refractive index of materials is elementary, yet it is crucial for the manipulation of 
electromagnetic waves. Relying on the inherent properties of natural materials, it has been a long-standing 
challenge in device engineering to increase the index-modulation contrast. Here, we demonstrate a 
significant amount of ultrafast index modulation by optically exciting non-equilibrium Dirac fermions in 
the graphene layer integrated onto a high-index metamaterial. Furthermore, an extremely-large electrical 

modulation of refractive index up to An 3.4 (at 0.69 THz) is achieved by electrical tuning of the density 

of the equilibrium Dirac fermion in the graphene metamaterial. This manifestation, otherwise remaining 
elusive in conventional semiconductor devices, fully exploits the characteristic ultrafast charge relaxation in 
graphene as well as the strong capacitive response of the metamaterial, both of which enable us to drastically 
increase the light-matter interaction of graphene and the corresponding index contrast in the graphene 
metamaterials. 

The control of refractive index in natural materials has attracted considerable attention in the field of 
holographic imaging systems, data storage, and for numerous photonic applications 1 . Changes of refractive 
index under external stimuli have been studied intensively in various transparent nonlinear crystals 2 ' 3 and 
light-sensitive materials 4 . However, these natural materials inherently have small index changes due to high-order 
nonlinearities known as Pockels or Kerr effects, or insignificant modification of molecular structures in polymers. 
More importantly, fast refractive index modulation in polymer-based materials also is limited significantly by the 
slow recovery time to the initial state. Therefore, promising alternatives to natural materials for ultrafast index 
modulation with large contrast are highly desirable. The emergence of artificially- engineered materials, referred 
to as metamaterials 5,6 , enables the manipulation of the electromagnetic properties of materials far beyond the 
limits of natural materials. By periodically arranging subwavelength-scale metallic elements, both negative and 
positive indices of refraction were demonstrated recently with extremely high values 7 ' 8 . More interestingly, the 
tuning ability in the electromagnetic responses in metamaterials, implanted by hybridization with various active 
media, recently has broadened the field of metamaterial applications 9 " 15 . 

Graphene, two-dimensionally arranged carbon atoms in honeycomb lattices, has unique optical absorption 
properties due to its linear dispersion and the massless nature of its carriers 16 . For example, a monolayer of 
graphene absorbs —2.3% of near infrared and visible light due to the universal optical conductivity of the massless 
Dirac fermion in graphene 17 . This constant and broadband optical absorption shows no Fermi level dependency 
except near the frequency threshold 18 for interband transitions ~2|£ F |. However, in the THz regime, where the 
optical conductivity of graphene is determined mainly by the intraband transitions 19 , the optical absorption 
spectra strongly depend on the Fermi level owing to its band structure. For this reason, the modulation of THz 
waves using a graphene layer recently was demonstrated successfully by electrically or optically controlling its 
Fermi level 19 " 21 . An atomically-thin graphene layer also has a well-known photoconductive nature 22 that can be 
explained by the diffusion of photo-excited hot carriers with a quasi-Fermi level 23 . The excited carriers undergo a 
relaxation process that is, in principle, ultrafast (—10 ps at room temperature) owing to the emission of hot 
optical phonons (phonon energy of w 0 = 196 meV), in which the energy loss is about an order of magnitude 
faster than that in conventional semiconductors, such as silicon (co 0 = 63 meV), germanium (co 0 = 37 meV), and 
gallium arsenide (co 0 — 36 meV) 24 " 27 . In addition to the ultrafast carrier dynamics of graphene, broadband 
operation is possible, because the real part of the optical conductivity has nearly non- dispersive values in the 
THz frequency range 19 . With this end in mind, we concluded that a graphene layer is a promising active medium 
when integrated with a metamaterial for ultrafast manipulation of broadband THz waves by electrically or 
optically tuning its refractive index. 
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Results 

Design of graphene metamaterials. To achieve active control of the 
refractive index in a THz metamaterial, an array of meta-atoms, an 
atomically thin graphene layer, and an array of metallic wire gate 
electrodes are functionally configured together into an ultra- 
compact, thin, and flexible polymeric substrate. The structure of 
the integrated THz graphene metamaterial designed for large 
index modulation is shown in Fig. la. As an atomically- thin active 
medium, a monolayer graphene layer was directly deposited onto an 
array of hexagonal metallic meta-atoms specially designed to have a 
high index of refraction in the THz frequency range 7 . As can be 
confirmed by the image of the saturated electric field in Fig. lc, 
large numbers of surface charges accumulate on the edge of the 
hexagonal metallic frame, which induces a significantly large 
polarization by strong capacitive coupling between the unit cells. 
Among the accumulated charges, some portion of the charge 
carriers leak into the attached graphene layer, resulting in a 
reduction of the induced polarization, which determines the initial 
refractive index of the graphene metamaterial at the charge neutral 
point (CNP). When the Fermi level of graphene is tuned by proper 
electrical gating, the induced polarization is altered because the 
intraband optical conductivity of graphene is highly sensitive to its 
Fermi level, and the conductivity change modifies the capacitance 
between the meta-atoms. With this principle of polarization change, 
the effective permittivity, and thus the refractive index of the 
graphene metamaterial can be electrically modulated by applying 
DC gate voltages. For the electrical tuning of the Fermi level of the 
graphene layer, a previously reported metallic wire array was used as 
a gate electrode 15 , which was devised to be transparent for the THz 



waves polarized normal to the wire direction. The transparent gate 
electrode also is beneficial because it has the possibility of applying a 
quasi-uniform DC electric field onto the graphene layer, in the same 
way as is traditionally done with thick and lossy silicon back-gates. 
An image of the fully-integrated device is shown in Figs, lb and d, 
and the geometrical parameters of the fabricated metamaterial are 
given in the caption of Fig. 1. 

Electrical modulation of refractive indices. Figures 2a and b show 
gate- controlled refractive index Re(n) and figure of merit FOM = 
Re(n)/lm(n) extracted from the THz-TDS measurement with an 
iterative algorithm considering multi-pass transmission 7 . At the 
CNP, a strong electrical resonance was observed with a peak refrac- 
tive index of 13.8 at a frequency of 0.58 THz and a refractive index of 
12.1 at the quasi-static limit where the loss of the metamaterial is 
minimal. We estimated the charge neutral gate voltage V C np of the 
device to be ~ 120 V from the measured trace of the transmitted THz 
field maxima (dashed line in Fig. 2). With increasing | AV|, where Ay 
— Vg — Vcnp an d Vg is the applied gate voltage, the LC resonance 
along with the peak refractive index were red- shifted and the 
resonance dip was diminished. The phenomenological resonance 
features of the graphene metamaterial arise from an increase in the 
real part of the optical conductivity in graphene. At a fixed frequency 
of 0.69 THz, the refractive index gradually decreased from 12.4 to 9.0 
with increasing |AV| due to a remarkable change of the capacitance 
between the meta-atoms. This unprecedented index contrast ARe(n) 

3.4 of the graphene metamaterial is much larger than that of any 

other natural material and hardly obtainable in nature. By 
conducting numerical simulation using the finite element method 




Figure 1 | Schematic view and optical micrograph image of high-index graphene metamaterials. (a) Schematic rendering of a graphene 
metamaterial shows that a monolayer graphene attached directly to hexagonal gold meta-atoms (a 60-um-size unit cell with a 3-um-wide, 100-nm-thick 
gold frame with a 5-um-wide gap between adjacent meta-atoms) and transparent wire electrodes (periodicity = 6 (am, metal width = 4 urn) are fully 
embedded in a thin and flexible polyimide (PI) substrate. The total thickness of the graphene metamaterial was —3.2 (am, which shows a high index of 
refraction (above 10) in the THz range. Both the ultrafast optical pump beam and electrical gating contributed to the modulation of the high index of 
refraction with large index contrast, (b) Optical image of the fabricated high-index graphene metamaterial. (c) In-plane electric field component 
(at 0.5 THz) of the graphene metamaterial at the CNP. (d) Optical micrograph of the fabricated metamaterial for a clear understanding of the realistic 
structure of the field image of (c). 
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Figure 2 | Electrical control of the refractive index and FOM of the graphene metamaterial. (a) Measured refractive index and (b) figure-of-merit (FOM 
= Re(n)/lm(n)) plotted as a function of gate voltage. With increased quasi-static carrier density, the refractive index decreases near the resonance 
frequency, (c) (d) Simulated results with a single-layer graphene approximation corresponding to a and b. Here, we approximated the carrier 
density at the conductivity minimum and intraband scattering time as 8 X 10 11 cm" 2 and 16 fs, respectively, (e) Gate-dependent refractive indices are 

presented at frequencies of 0.69 THz and 0.58 THz. At these two frequencies, a maximum index of 14.3 and a maximum index contrast of 3.4 are 

clearly shown, (f) Frequency- dependent refractive index profile with respect to the two different gate voltages, V cnp and V ma ^. The numbers in parentheses 
indicate the values of gate voltage. 



(FEM), the refractive index and FOM of the proposed metamaterial 
were reproduced and matched well with the experimental results, as 
shown in Figs. 2c and d. Considering carrier density at the conduc- 
tivity minimum 15 ' 28 in CVD-grown graphene to be n min = 8 X 
10 11 cm" 2 , the measured refractive index was in excellent agree- 
ment with the simulated values over the broadband frequencies. 

Optical modulation of refractive indices. To investigate optically- 
controlled ultrafast index modulation in the graphene metamaterial, 
we used ultrafast optical-pump THz-probe spectroscopy (See 
Methods and Supplementary Information). As described in the 
previous paragraph, the capability of Fermi-level tuning by 
electrical gating is quite unique, because it enables us to perform 
ultrafast time-resolved index modulation under various quasi- 
static operating conditions. As shown in Fig. 3, we measured the 
pump-induced THz electric field changes AE(t) for each fixed 
pump-probe delay At along with optically unexcited THz field 
E 0 (t). Fourier transform of the measured fields has the information 



on both real and imaginary values of the pump-induced As(co) or 
Acr(co). Here, As(co) or Ag(co) denotes the change of dielectric 
constant or conductivity in the active layer, i.e., meta-atom/ 
graphene layer (For the meta- atom-free samples shown in the 
Supplementary Information, the values are the change for the 
graphene layer). As displayed in Fig. 3a, the pump-induced AE(t) 
traces showed substantial field reshaping compared to the E 0 (t). First, 
we observed negligible phase shifts of AE(t) in the low-frequency part 
of the early THz field delay (t < 0), and apparent phase shifts in the 
high-frequency component of the THz field (t > 0). Second, the 
opposite sign of AE(t) is significant when t > 0, but it is scarcely 
appreciable when t < 0. As discussed later in more detail, these 
examinations indicate a strong in -phase resonant response at low 
THz frequencies and an inductive Drude-like behaviour in the high 
THz frequency range. 

The resonant in-phase response was characterized further to cor- 
roborate the electrical index modulation in the previous section. As 
shown in Fig. 3b, when the pump fluence was 18.7 |aj/cm 2 , for 
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Figure 3 | Ultrafast optical modulation of the refractive index in the graphene metamaterial. (a) Schematic of optical-pump THz-probe spectroscopy. 
The THz field through the graphene metamaterial without a pump E 0 (t) (gray dot) and with a pump (black line) are shown at different pump-probe 
delays At (b) Gate-dependent refractive index changes at At = 1 ps. The large index contrast ARe( n) — —2.4 and the red-shift of the resonance are clearly 
visible, (c) Photo-excited conductivity changes (red circle) in the active graphene layer of the metamaterial at a pump fluence of 18.7 uj/cm 2 near the 
CNP. The graphene layer shows a negative sign of the differential conductivity change near the resonance frequency, (d) Conductivity changes at 
0.83 THz and 1.8 THz show opposite sign, (e) Resonance frequency (parameters from Drude-Lorentz model) was slightly blue- shifted from At = 1 to 
4 ps due to the decrease of the graphene's intraband optical conductivity. 



example, the photo-induced change in the refractive index showed a 
large index contrast ARe(n) — —2.4 at around 0.61 THz, and the 
modulation depth was highly sensitive to the gate voltages that were 
applied. We note that this large negative-index modulation is hardly 
observable in any natural photorefractive materials or nonlinear 
crystals. In addition, the red-shift of peak refractive index with an 
applied gate voltage was clearly resolved as the gated Fermi level was 
gradually detuned from the CNP. This behaviour is due to the 
increased THz conductivity, and it is consistent with the previously 
described electrical index modulation. 

To substantiate our observations, we present the transient 
dynamics of Act (cx>) for each pump delay At at the CNP in Fig. 3c. 
Immediately after photo -excitation, the real part of conductivity 
Acr^co) exhibited an asymmetric frequency- dependent reduction 
Aai(co) < 0 near resonance frequency. Likewise, the imaginary part 
Act 2 (co) showed complementary zero-crossing, exactly at the res- 
onance frequency. The signs and shapes of these conductivity spectra 
indicate an increased on-resonance transmission due to the increase 
of leakage current into the graphene layer and the corresponding 
reduction of the capacitive coupling between the meta-atoms 
(Fig. 3d). These features were further confirmed by examining the 
time-dependent shifts of the resonance frequency as a function of 
pump-probe delays (Fig. 3e). During the carrier relaxation process 
(At = 1 ~ 4 ps), the resonance was slightly shifted to higher 



frequency because the carrier recombination process induces Acr^co) 
to decrease, and consequently leads to an increased capacitive- 
coupling between the adjacent meta-atoms. This behaviour agrees 
qualitatively with the time- domain in-phase analysis described above 
as well as with the electrically-controlled index changes in the previous 
section. In addition, this distinct resonant behaviour was ensured by 
performing the same experiment and checking Act (co) on the same 
sample without the metamaterial as a control sample (See Supple- 
mentary Information for the change in the non- equilibrium conduc- 
tivity of the control sample). As expected, no strong resonance was 
observed over the entire frequency range; the measured spectra were 
featureless with broadened Act^co) by the grain-boundary-limited 
charge back-scattering in the Drude response 29 . 

Ultrafast carrier dynamics of graphene metamaterials. To under- 
stand the nature of ultrafast index modulation, we investigated the 
gate- dependent non-equilibrium THz dynamics. Figure 4a is a plot 
of the negative transmission change —AT(t)/T = —[2AE(t)/E 0 + 
(AE(t)/E 0 ) 2 ] at a fixed field delay at t = 0 as a function of the applied 
gate-voltage. As discussed above, the in-phase response (no phase 
shift) at early field delay justifies that the overall dynamics is 
determined by measuring the changes in peak AE(t). The maxi- 
mum AE(t) occured within the time-resolution of our THz pulse 
(~1 ps), and the relaxation signal contained multi- exponential 
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Figure 4 | Ultrafast dynamics of photoexcited carriers in the graphene metamaterial. (a) Negative transmission change ( — A 77 T 0 ) at 1 8.7 uj/cm 2 pump 
fluence at three different gate voltages (A V = 0 V, — 100 V, —200 V). The measured values are fitted by a solid line using a bi-exponential decay 
function. Inset: The band diagram of a monolayer graphene with gate-dependent Fermi levels, (b) Enlarged view of (a) with a logarithmic scale of the 
transmission. The two fitting lines at each gate voltage intersected at different pump-probe delays. We attributed the slow decay to the acoustic phonon 
scattering. Early emergence of the slow decay with increasing gate voltage is clearly shown. 



decaying components. Note that the data showed a strong gate- 
dependent feature for the peak AE(t) changes. In contrast to the 
interband optical spectroscopy, where the transmission changes of 
the optical probe pulse (typically >2|£ F |) are dominated largely by 
the photo-excited carrier density, the THz spectroscopy is highly 
sensitive to the initial Fermi level of graphene, since it read the 
carrier distribution in all energy spectra, including the initial 
carrier density as well as the photo -excited carrier density. As 
shown in Fig. 4b, more detailed analysis exhibited that a bi- 
exponential fit using (1 — B)exp( — t/Ti) + (B)exp( — £/t 2 ) properly 
described the THz transient dynamics; we found that a mono- 
exponential fit did not reproduce the data. Using a logarithmic 
scale, the data seemingly represented two different time scales. We 
attributed the picosecond fast relaxation component x x to the optical 
phonon emission (co 0 = 196 meV) of the excited hot Dirac fermions 
and the other t 2 to the slow acoustic phonon scattering of low-energy 
carriers 23 ' 26,27 ' 30 . Combining the two relaxation components with the 
gate- dependent dynamics, we observed that both the amplitude B 
and the corresponding scattering rate 1/t 2 of the slow relaxation rate 
became larger with increasing |AV|. This implied that the effect of 
acoustic phonon scattering emerged earlier under high carrier 
density 26,31 . Given the photo-excited carrier density (—1.2 X 
10 12 cm" 2 ), it is well known that the initially photo-excited hot 
Dirac fermions share their energy with the existing cold carriers 
(densities determined by |AV|) through an ultrafast thermalisation 
process. As |AV| increased, the initial carrier temperature after 
thermalisation decreased so that the excited carriers had less 
chance to emit multiple optical phonons 26 , thereby exhibiting the 
observed early emergence of acoustic phonon cooling as the 
dominant relaxation mechanism. Thus, the ultimate speed of 
ultrafast index modulation was limited largely by these persistent 
cooling dynamics of hot Dirac fermions. 

Discussion 

By the electrical and/or optical tuning of the Fermi level of monolayer 
graphene, capacitive coupling between unit cells of graphene- 
attached metamaterial can be controlled, causing a corresponding 
change of the effective index of refraction. Based on this principle, we 
designed high-index THz metamaterials integrated with a CVD- 
grown monolayer graphene, and transparent electrodes embedded 
in an optically thin (—2/150) flexible polyimide substrate. First, 
electrical modulation showed that an extremely large index contrast 



of ARe(rc) 3.4 was achievable simply by applying a gate voltage, 

and the observed value was unprecedented and hardly found in 
nature and can be increased further by using more advanced gating 
techniques. Second, we have provided compelling evidence that sup- 
ports the possibility of achieving ultrafast optical control of the 
refractive index at a picosecond time scale. Spectrally- and tempor- 
ally-resolved THz dynamics show that the strong coupling between 
the optically- excited graphene and the resonance of the metama- 
terial led to optically- controlled capacitive coupling in the graphene 
metamaterial, manifested by an increase of the on-resonance trans- 
mission. Third, we have presented a detailed investigation of the 
gate- dependent ultrafast Dirac fermion relaxation dynamics. By elec- 
trically changing the Fermi-level, we can understand that optical 
phonon scattering and acoustic phonon scattering are the key lim- 
iting factors in achieving ultra-high- speed index modulation. 
However, we emphasize that the relaxation process was significantly 
faster than that of almost any other semiconductor devices. Thus, 
graphene, if used as an active layer, may have a great advantage in the 
ultrafast operation of various THz devices, as well as in the index 
modulation of THz metamaterials for tunable transformation optics. 

Methods 

Fabrication of the graphene metamaterial. All metallic parts of the refractive index 
tunable metamaterial were made of 100-nm-thick gold with a 10-nm-thick 
chromium adhesion layer. The construction of the grapheme metamaterial was 
initiated by the spin-coating of a polyimide solution (PI-2610, HD MicroSystems) 
onto a sacrificial silicon wafer. After a two-step baking process, the polyimide film was 
fully cured, resulting in a final thickness of 1 urn. To define the transparent wire 
electrode, UV photolithography and electron-beam evaporation were followed by a 
metal lift-off technique. The metallic meta-atoms were patterned by repeating the 
same polyimide stacking processes and metal patterning as employed for the first 
layer. Then, commercially- available graphene from Graphene Square, Inc., grown by 
CVD onto copper film, was transferred to the entire area of the metamaterial (10 X 
10 mm 2 ). After the graphene transfer, a ground electrode was defined on the 
graphene layer using a shadow mask, and we finally covered the device with a 1-um- 
thick polyimide film. After opening the electrical contact via oxygen plasma etching of 
the covered polyimide film, the thin, substrate-free, flexible tunable graphene 
metamaterials were peeled off of the silicon substrate, as shown in Fig. lb. Finally, the 
active graphene metamaterial was soldered to a drilled PCB substrate. 

Optical-pump THz-probe spectroscopy. To characterize the electrical and optical 
modulation of the refractive index of the fabricated graphene metamaterial, we 
conducted optical-pump THz-probe spectroscopy together with electrical gating. The 
graphene metamaterial was excited by a nearly bandwidth-limited ultrashort 50 fs 
pulse with a photon energy of 1.55 eV delivered from a 250-kHz Ti:sapphire 
regenerative amplifier (Coherent RegA 9050). A fraction of the amplifier output was 
used to generate the THz pulse via optical rectification in a 500-um-thick <110> 
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ZnTe nonlinear crystal. The collimated THz field was focused with a two-inch gold- 
coated off-axis parabolic mirror onto the high index graphene metamaterial. To 
detect the transmitted THz field through the sample, we used field-resolved electro- 
optic sampling in the same ZnTe crystal. The measured transient THz response 
directly reflected the dynamics of the complex Ae(<jj) or Aa(co), and the two are 
related by Ae(co) = As^eo) + iAs 2 {co) = — Ao" 2 (co)/co£ 0 + /A<Ti(co)/coe 0 . All 
measurements were performed at room temperature. 
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